ON THE DISTRIBUTION OF THE FREE PATH LENGTH OF 
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Abstract. Let C > 2 be an integer. For each s e (0, i) remove from R 2 the union of discs of radius e 
centered at the integer lattice points (m, n), with m £ n (mod ('). Consider a point-like particle moving 
linearly at unit speed, with velocity ai, along a trajectory starting at the origin, and its free path length 
Ty E (w) e [0, oo], We prove the weak convergence of the probability measures associated with the random 
variables etc .e as e — > + and explicitly compute the limiting distribution. For I = 3 this leads to an 
asymptotic formula for the length of the trajectory of a billiard in a regular hexagon, starting at the center, 
with circular pockets of radius e — » + removed from the corners. For t = 2 this corresponds to the 
trajectory of a billiard in a unit square with circular pockets removed from the corners and trajectory starting 
at the center of the square. The limiting probability measures on [0, oo) have a tail at infinity, which contrasts 
with the case of a square with pockets and trajectory starting from one of the corners, where the limiting 
probability measure has compact support. 



1. Introduction 

Recent progress led to a better understanding of the statistics of the free path length of the periodic 
Lorentz gas in the small scatterer limit [2,3,5,6,7,8,9,10). The aim of this paper is to study situations 
where periodicity conditions are altered by imposing certain congruence conditions on the integer lattice 
points where scatterers are placed. The case of the honeycomb lattice arises as a particular example in 
this context. Here we only consider the situation where motion originates at the origin, extending the 
results from |2] and [3|. The case where the initial position is randomly chosen is more intricate and 
will not be treated here. 

Let i > 2 be an integer. Consider the set Z 2 ^ of pairs of integers (m, n) with m £ n (mod £), For 
every e > consider the "fat lattice points" (scatterers) given by small discs of radius s centered at all 



points of and the region 

Z ts = {x e R 2 : dist(x, Z 2 Q ) > s] 

obtained by removing all scatterers. In R 2 consider a point-like particle moving at constant unit speed 
along a linear trajectory originating at (0, 0). The free path length (first exit time) is defined as 

T( ye (u>) = inf{r > : toj e dZ^ s ), 

the distance traveled to reach the first scatterer along the direction to = e"° e T, to e [0, 2n), and as 
+co when the particle escapes to infinity without reaching any scatterer. The Lebesgue measure of a 
measurable set A c R is denoted by |A|. This paper is concerned with the study, in the small scatterer 
limit (e — > + ), of the asymptotic behavior of the repartition function of ST( >S defined by 

{we[0,27r):T f , £ (w)>^} 



w - 1 



To accomplish this we first consider the situation where scatterers are obtained by translating the 

■2 
V) 



vertical segment V e = {0} X [s,s] by (m,n) e ZjL and estimate, for any interval / c [0, 1] of length 



x s c with fixed c e (0, 1), the repartition 



G e ,i,eU) 



cd e arctan I : <jy, e («) > - 



+ 
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of the horizontal free path length 

q{, £ (<jL>) = inf {g : (q, q tan at) e Z 2 ^ + V £ J , we 



0,: 



In this paper tp will denote Euler's totient function. The dilogarithm is defined by 

ln(l -t) 



Li 2( x) = £ x ~ 2 = - r . 



■ off, 



xe [0,1]. 



Clearly Li 2 (l) = f (2) = £. 

The main result of this paper shows that the limit of P^ >e exists as s — > + and this limit is explicitly 
computed. 



Theorem 1. (i) For every < cj < 1 and 5 > 0, as s — > + , 

G f ,/, e (/l) = C/ G^i) + ( e - 5+fl ( c - c >') , A > 0, 

w/iere 

c ' - J I d + u 2 ' 6>(c,ci) = min|c + ci, i -2ci|, 
one/ f/ie limiting repartition function G[ is given by 

■\-(JL-. + A{t))A ifAe(0,i], 



l-^+AW// 2 (l) if^6[i,l], 



2C(0 



if ie [l,oo), 



with 



1 2C(f) 



/? prime 



p prime 



H 2 (A) = 3A - 2 + ({2) - (In Af + 2(1 - A) In | - - 1 1 - 2 Li 2 (^), 
ff 3 (A) = LiJi)-(A-l)lJl-~)-l. 



(ii) For every 6 > 0, as s — > + 



f , e (^) = GKi) + O w (^+' s ), ^>0. 



(1.1) 



(1.2) 





Figure 1. 77ie repartition function G3 anrf ?/ie density function g$ 
The continuity of Gf at A — I is equivalent with the well known dilogarithm identity 



A 1 /n r 1/2 i 1 



<T(2)-(ln2) 2 
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Since ^ ^ and A(() -> jL: as t — > oo, the compactly supported limiting repartition //(/I) from 
|H] Theorem 1.1] is being recovered as lim^oo G((A). 




Figure 2. The free path in a hexagonal billiard and in a honeycomb 

Our original motivation for considering this problem comes from the study of the exit time of the 
linear motion with specular cushion collisions on a hexagonal (open) billiard table with (small) circular 
open pockets of radius s removed from its corners (see Figure©. The starting remark here is that, after 
unfolding the hexagon to a honeycomb in R 2 , one can deform the later to Z 2 3) . This process converts 
the problem on the hexagonal billiard with pockets into one concerning the free path length of a Lorentz 
gas in R 2 with small identical ellipses centered at the points from Z 2 3) as scatterers (see Figure[8). 

Let Tg (to) denote the free path length in the hexagonal billiard with discs of radius s removed from 
the corners and motion starting at the center, and let 

A 



Le [0,2tt] : r f hex (w) > - 

denote the repartition function of er] lex . We prove 
Theorem 2. For every 6 > 0, as e — > + , 

P e hex (A) = G 3 + O s (e*~ s ) , A > 0. 

Scaling e to e V2 one can apply Theorem [T] (ii) with t = 2 to estimate the repartition function 
¥°(A) = [0, 2n] : t°(w) > ^}| of the free path length t°(oj) of a billiard in the unit square with 

pockets of radius s at the corners and trajectory starting at the center (see FigureQ, getting 

Theorem 3. For every 6 > 0, as s — » + , 

V n JX) = G 2 (^ + O s (s L *- s ), A>0. 

Theorem [3] should be compared with the situation where the initial position is at one of the four 
vertices, and where the limiting distribution has compact support (2] |3]. It is not clear whether these 
methods would directly extend to other concrete initial positions, such as (i,0), n e N. However, it 
looks likely that further refinements could lead to "space-phase average" results similar to those from 
131 . The main difficulty seems to arise from the increasing complexity and number of cases that need 
to be analyzed in detail, leading to integrals in the main terms of the asymptotic formula which are 
manifestly more intricate than in the case of the square. 
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Figure 3 . The free path in a square billiard with trajectory starting at the center 



2. The contribution of consecutive Farey fractions y < y' with y, y' e Tq \ T^ 

The integer part of a real number x is denoted by [x] . Let s > and Q = [-]. Denote by Tq the set of 
Farey fractions y = - in lowest terms with < a < q < Q. For any interval / c [0, 1], set Ti,q = IC\Tq. 



Denote by Tq the set of Farey fractions y — - e Tq with ( \ (q - a), and set Tjq -In T Q . Set also 
jr(C) _ Ug =1 !Fg^. It is well known that if y = - < y* = 4 are consecutive elements in Tq if and only if 

a'q - aq' 

In particular we have 



1 and q + q > Q > max{g, g'} 
e(<7 + g') > e(g + !)>!> max{e^, eg'}. 



(2.1) 



For convenience consider 

Gt.iJA) := \{co G arctan/ : qe, B (co) > AQ}\ 



[As 



A>0. 



For any interval / c [0, 1] with |/| x Q e , < c < 1, we will prove a formula of type 

G^ S (A) = c 7 G,(i) + 5 ( e -Wi>) , 

which will immediately imply dl.lt . because e[i] = 1 + 0(e) and for any compact set K c 
C# > such that 

\Gdx') - G t {x")\ < CkW ~ 
Given y < y' consecutive elements in Tq, denote to = y 
(|2~TT i and 

1 - sq s(q + q') - 1 , 1 - eq' 



there is 

x',x"eK. (2.2) 
f = 2^£, M = r + f = Employing 



?o - y 



93 



MO - to 



qq 



y - m 



33 



(2.3) 



we find y < ?o < mo < y'. In particular if y < y' < y" are consecutive elements in Tq, then y + | < y' < 
y" — %. Therefore the intervals [y--,y+-],y= 2 £ cover the interval [0, 1] in such a way that 
every element in [0, 1] belongs to at most two of these intervals. As a result any trajectory with slope 
tanoi e (y, y') will intersect, as in the case of the square lattice flUE], one of the scatterers (q,a) + V e or 
{q',a') + V s . 

Since a'q - aq' = 1, only two situations can occur here: y i T^ and y' £ T (t> , respectively y e T {t) 
or y' e !F (f) . The later will be discussed in Section 3. In the remainder of this section we assume that 
y £ T (C) and y' i. T^\ situation where the horizontal free path is given (see FigureE} by 



3 

min{<7, q'} 

q' 



if y < tanw < to, 
if to < tanw < mo, 
if mo < tan co < y' . 
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qcjw) = q 



qcju)) = mh%, q'} qt, E (aS) = q' 

V 



r t uo y 

Figure 4. The horizontal free path when y, y' e Ti,q \ 



In this way the contribution of the interval [y, y'] to Ge,i,e(A) is given by 

if max{q,q'} < AQ, 

arctany' - arctany if min{q,q'} > AQ, 

arctany' - arctan«o if q < AQ < q', 

arctanfo - arctany if q' < AQ < q. 

This contribution is zero whenever A > 1, so we next assume < A < 1. 
Using (12.3b . the estimates 

h , 1 1 , 2 

arctan(x + /i)-arctanjt= - + 0(/r), r - < 2(y' - y) = — , (2.4) 

1 + x 2 1 +y 2 1 + y' 2 qq' 



and the inequality 



Z « 2 «' 2 ^ Z 



■ y i + y- 
1 1 



we infer that the contribution to G^/^Ol) of all intervals [y, y'] c / with y < y' consecutive elements in 
Tq and y, y' £ IF^ is given by 



with 



}fl s (A) = A /ig (A) + B m (A) + C 7 , e (i) + O(s), 

1 1 



A/, fi (A) = V — ■ - — 
^ Flf) qq' i + r 2 



min(^,^'|>/12 



y,y'eT hQ \r 
q<AQ<q' 



l-sq' 1_ 

qq' i+r' 



q'<AQ<q 



1 — eq 1 
1+y 2 



Since there are no consecutive elements in Tq which belong both to T { ^ we have 



Ai,q{X) = Al Q (A) - AJ£(A) - AJ*(X) 



(2.5) 



(2.6) 



(2.7) 



with 



1 1 



rri ff<7 1 



min{iy,g'J>/l<2 



t/g 1+y 



2 ' 



Z 

min{#,#' |>/12 



l l 



1+y 2 ' 



A /i 2 ^ = z 

min{#,#'}>,l<2 



l l 



re r,, ,/<> qq 



1+y 2 ' 



and respectively 



= B + IQ {A) - B-/JX) - B-fl{A), C,, Q (A) = C+(A) - Cf {A) - CJ^(A), (2.8) 
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with 

q<AQ<q' q' KAQ<q 

**o= Z ^-A- c ^ )= Z ^tA. 

e ,^n a 1+7 " 2 ,<£W e qq 1+7 

q<,AQ<q' q'^AQ<q 

Lemma 1. For any function V e C'[0, iV] w/f/i fofaZ variation T^V and C(i) as in fll.2| >, 

V ^V(t)=C(0 f V + O c ((\\V\U + T^V)lnN). 
« Jo v 



g cd(a)=i 



Proof. The left-hand side L^jy can be expressed as 

Z Z^« 

g cd(a>=i 

Writing k = dfc' with 1 < A 7 < ] and gcd(^of) = gcd(^fe') = 1, denoting V n (x) = V(nx), (Tq{€) := 
#{d > 1 : d | £}, and using Mobius summation (as in Q] Lemma 2.2]) we have 

^- Z f Z wo 

gcd(f,rf)=l gcd(f,/t')=l 

-[N/d] 

h \ r I, 

gcd(W)=l 



gcd(W)=l 

Z ^-3(rv + wviu)) + 4iviu + 7yv)Zl 

1 l«d<2V « " WO / \ d=l " 

gcd(M)=l 

y ^ ^ f V + O(||V|UlniV) + O^(0lV|U+7^V)]n^) 

gcd(f,d)=l 

z ^^(z^))r y+ ^ (||y|u+r ' y)in ^ 



d»l d>N 
gcd(f,rf)=l 



with 



= C(i) V + 0,((||V|L+r*V)liitf), 

c<) " « & « 2 ~ ' LH W'iw'Jl 1 

p pnme p pnme 

which gives the desired estimate. 

Lemma 2. For any function V £ C l [0,N] and any S > 0, 

Y j n-lv(n) = €C(€) J o V + O e , s ((\\V\L + T^V)N s ). 
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Proof. Let t = pV ■ ■ ■ p" r with p\ , . . . , p r distinct primes and ai, . . . , a r e N (so r = co({), the number 



'i 

of prime divisors of €). Writing n = p 1 ^ ■ ■ ■ p kr m with k, > and gcd((, m) = 1, we obtain 

h+a:- 



v y( fa ) T/ / n v n ^ 2 v t/ / \ , ™ 

y(n) = z, 1.1—?— ^ ^r^-^ (m) - (Z9) 



l<iXiV fa,..,*,*) l<S<r Pi Um ^ /( ^i...^) 

Pi-Pr r <N gcd(f,m)=l 



According to Lemma[T]the inner sum above can be expressed as 

r Ni(p\ i -pin 

C{€) J Q V{p k ; ■ ■ ■ p k 'x) dx + O e {{\\V\U + T%V) lniv) 



- k C{{) 1 f V + < ((||V|U + rS r V) In at). 
Pj 1 ■ Jo 



Inserting this into (|2.9l l and using the fact that the number of terms in the first sum in (12. 9t is <k (In N) wi - e) 
and (p(p ki+a ') = p ki+a '(l - pj 1 ) we infer that the expression in (|2.9b is given by 

^n(i--l y - 1 , f y + o f ((||y|u + ^y)(in^v) 1+ttOT ). (2.10) 

p prime p k , 1 -pf <N 



The statement now follows from ( 12.1 OK using 



fa,..,*,>0 ^1 Pr p|f v ' 7 

/? prime 



and the bound 



I 



(In AO 



r-l 



, , fa A-,. "•PU-.Pr N 

h~>0 Pl '-'Pr " 

a 

The following estimate JU Proposition A4] will be employed several times. 

Lemma 3. Assume that q > 1 and h are two given integers, I and Sf are intervals of length less than q, 
and f : I X Sf — * R is a C 1 function. Then for any integer T > 1 and any 5 > 

y f{a,b)=^ ff f(x,y)dxdy + 8, 

(mod 
gcd(&,#)= 1 



with 

G« s r-'n/iL.^'-gccK/,.,/)- -i niv/iL,./- ■'»,-d(/j. f/ )T + 

where llflloo l|V 



1^1 + 


1*9/1 


1 3a- 1 ^ 


1 3v 1 



Lemma|3]will be typically applied to the following situations: Let / be a subinterval of [0, 1]. For 
every q e [l,Q] consider the intervals I — ql and J' - J x, q - (max{AQ, Q-q},Q], and the functions f q 
and g q defined on ql x J x, q by 

11 1 - sv 1 

f q (u,v):= gq(u,v)= —z. 

qv 1 + qv 1 + 

We clearly have 

H/slL^i HV/,|U «, Woc«^, HV^IU «, ^-L (2.11) 
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Proposition 1. For every A e (0, 1] and c\ e (0, 1), 

where 



f'i 

Ja x 



maxRl-*} ((hU) 2 j/A€[i,l], 
and 9\(ci) = max {2ci - j, -c\, -lj. 

Proof. There is at most one y e fi,Q with y' i. I. Since j|| «: Q c and < the total contribution 
to the final asymptotic results from Theorem[T]resulting from replacing the two conditions y, y' e Ti,q 
by a € ql or by a' € q'l will be <K (2 1 and respectively <sc Qs -1 , thus negligible. As a result we shall 
tacitly do this in formulas (12121 . ( 12.151 . (12191 . (I%24l and (l2~27l 

Furthermore, the summation constraints in 7 = 2 € T^q and / = ^ e !Tg will translate, using 
standard properties of Farey fractions, into the following constraints on the triplet (q, a, q')\ 

iq e {AQ, Q], 

\a€ql, q' € J^ q , gcd(q' , q) — 1, aq' — -1 (mod g). 

Summing first over q and then after integer pairs (a, q') e ql x J^ q as above (note that the number of 
such pairs is < q\I\) we infer 

A i q W= Z Z + (2.12) 

gcd(<j',f/)=l 
«^'=-l (mod q) 

Applying Lemma[3]with T = [Q C[ ] and (12. lit , the inner sum in (12.12b can be expressed as 
<7 2 Jy,, «v J ?/ : + ^) 2 v y 

ln — ,m ?i + °w (e 2ci ^ +5 + e _1_ci ) ■ 



(2.13) 



9(4) , 1 

= C[- 

q~ max (/1, 1 - Q , 

Summing over q e (/l(2, £>] in ( 12.13b and employing (12.12b and [ 1 , Lemma 2.3] we find 

Al Q {A) = ^ h{A) + O s ,a (Q s+e ^) . (2.14) 

To estimate Afg(A), note first that gcd({, q) — 1 because t\q — a and gcd(g, a) — I. As a result, 
putting w — q - a = (u, v — q' and using a'q - aq' = 1 and q' > Q - g, we find (with I denoting the 
multiplicative inverse of I (mod q)) 

Z Z ±-nsH* 

(2.15) 



^G<9<G (H>,v)e 9 (l-/)x7,,., 1 + (2-^) 
gcd(f,9)=l flw, gcd(v,?)=l 
wv=l (mod </) 

- Z Z /,(<Z-^v) + o(i). 

gcd(f, 9 )=l gcd(v,?)=l 
wv-^ (mod 

By Lemma|3](with T = [Q C[ ]) and ( 12.11b the inner sum in ( 12.151 ) can be expressed as 

'"- max [A, 1 - o j 



(2.16) 
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Summing over q in ( 12.16b and employing (12.151 and Lemma[T]we find 

A^H) = Cj ^UX) + O SM ((?^). (2.17) 
To estimate A~' 2 we use the second inequality in ( |2.4t and Syey, e ^ ^ 1 t0 infer 

A ?n= y — • 1 , + 0\-). (2.18) 

mm{q,q']>AQ 

We now proceed as for AJg, noting that gcd(^, g') = 1 because ^ | q' -a' and gcd(g', a') = 1, and getting 
as above (with w = q' — a' = ^m, v = g) 

a;^)= y y -i 1 +0 (± 

?'e(iQ,G] (»M0e 9 '(l-/)x7 v * 1+1— ^ I ^ 
gcd(^')=l f| w ,gcd(v, 9 ')=l 
vvv=-l (mod 



Z Z w-^+ofi) 



gcd(f, 9 ')=l gcd(v, 9 ')=l 

uv=— { (mod g') 

qe{AQ,Q] \ * UJ *A' * ° f (1 'J 1 + l q ' J 
gcd(f, 9 ')=l 

££ y ^) ln 5 — ^ + o s[ (q^-^ + q--) 



q'e(AQ,Q] * < ) I 

gcd(f, 9 ')=l 



c/C(0 



/iW + o w (e 5+fli(ci) ). 



The estimate for A />e (i) follows from < |2T7l d2~14l i, ( f2~T7] > and d2~18l i. 
Proposition 2. For every /I e (0, 1] and c\ e (0, 1), 

B I>Q (A) = c,A(%(A) + Oftw (Q s+S ^) and C m {X) = C/ A(€)/ 2 (A) + O SM {q 5+( >^) , 

where 

r 1 i-x a 

hU) = In dx. 

Jmax{A,\-A\ x I - X 

Proof. Using the argument leading to ( 12.181 ) we see that 

q<AQ<q' 

Using customary properties of Farey fractions we infer (setting u — q - a, v — q') 

q<AQ(u,v)eq(l-I)XJ Xl , 
uv=l (mod q) 

Applying Lemma[3]to g q and T - [Q C] ] to the inner sum above and using (12.1 It we find 

bm= y W f f l -^d V+ oj Q ^ h+s ■ QCi ^ s ■ q2 ^ A 



q<AQ q 
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where 

V(n) = - f — -dy, ne[l,AQ]. 

n JmaxRl-g) y 

Using 



<± and rfy^i (2.20) 



and applying Mobius summation to V (e.g. 0] Lemma 2.3]) and Tonelli's theorem we find 

d f AQ du f 1 1 -y 

'max{A,l- 



M £(2) Jo U Jnux/j )> V ' 



e 



=t^ r — r ^^+o*(c2 wi<w ) (z2i) 

f(2) Jo * JmaxRl-x} J V 7 

= ^/ 2 w + ^(fi a+fli(ci) ). 

The condition £ \ q - a gives gcd(^, q) — 1 . Taking q - a - (u and v = g' and proceeding as in the 
case of Aj'q from the proof of Proposition [T] we have 

B m^= Z ^'177 

= ( U , v o e(? /|5-/)xy.,^'l + (^) 2+0 (e) (2.22) 

gcd(£ 9 )=l gcd(v,?)=l y 
«V=t (mod q) 

c i v ^(?) 



gcd(£,«)=l 

Applying Lemma[T]to the last sum in ( 12.221 1 we find 

ciC(€) C AQ du C 1 l-y 



JO u JmaxRl-*) y 



C,C{€) 



y (2.23) 



To estimate B l '1 we first fix q e (AQ, Q], then set u — q-a, v — q'. Since l\q' -a! and gcd(g', a') = 1 
we have gcd(^, q') — 1 . Moreover, q' - a' — tZZ— is divisible by I, somve 1 (mod (fq)) and we have 

i_ q 

B ^ = Z ^'^^Z Z ft(«-".v) + 0(i). (2.24) 

gs^U^ (mo(J ^ 

By d2.24| >, Lemma|3]and d2.1 It we infer 

(<p(i q ) 1 r 1-5 



= g^^) v(9) + 05 ( GI nax (2ci -H,- Cl) ) 

' q<AQ 1 



(2.25) 



From ( 12.25b . Lemma|2]and (12.20b we infer 



fljJO*) - ^ X V + W (e**«») = ^ h(A) + S/ (e' 5+fl ' (Cl) ) . (2.26) 
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The desired estimate on Bi,q(A) follows from d2.81 >. ( 12.21b . ( 12.23b and ( 12.251 . 
In similar fashion one gets 

c ^ = Z ^ - rr^ = Z Z ^ - «. v) + o (i) 

q'<:AQ<q %cd{v,q')=\, uv=-\ (mod 4') 

= ^ / 2 M) + o J (e i *'"" 1 ). 
<3»- Z Z «rtf-^ + o(i 

yeF, (f > r'en e 9 ' </lG (".'')=(?' -«',9)e?'(l-/)x^y xii 

q'<AQ<q gcd(v,?')=l KV=-1 (mod (V)) 



(2.27) 



c 3«- Z ^ ttW= 2 Z *w-*« + o(i) 

q'<AQ<q g"V,«J gc( J(v, ? ')=l,BV=-^(mod ? ') 



-jmax{2ci - ^ +(S,-ci — 1 
gcd(f,«')=l 

= ^W) + w (<2 tf+Wci) ). 

The desired estimate on C/,g(/l) follows from ( 12.8b and ( 12.27b . 
Corollary 1. For every A > ant/ 5 > 0, 

G« e (A) = C/ A(^)G (1) (i) + Oy( £ - s+8fec '»), 

where 

ln(l - A) In A + i In £ d« + 2 JjV, ±=s In £ <i« (0, ±], 

M 2 + 2jfl=sin T i ; Jii «T A e [|, 1], 

if A e [I, oo). 



G m (A) 



3. The contribution of consecutive Farey fractions y < y' with y e !Fg f) or y' e Tq 
Suppose first that y < y' are consecutive in Ti,q and y e where again Q — [i]. Consider 

ak = ka + a , qu = kq + q , y^ = — , tk — yk = , k > 0, t-\ = y , 

qk qk qk 

Inequalities ( 12.1b show that 

k 

y < — yk+\ < Jk < ■■■ < 7i < 7o = V and t k+ \ <t k < y k +i < yk- (3.1) 

Since q = a (mod €) and a'q — aq' — 1, we cannot have q' = a' (mod €), and so y' e Ti,q \ ■ 
Moreover, since I \ (q - a) we must have gcd(^, q) — 1 and qk £ ak (mod €) for all k > 0. 



The number K — [ ' — q ] is the unique integer > 1 for which q^ = Kq + q' < /IQ < qic 
(K + l)q + q' . The presence of the sink at y, (13.11 ) and yk t show that (see also Figure|5]) 

9f,e( w ) = 1k+i if f*+i < tan w < t k , k > -1, 



+ i 
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qiA^) = 1k+i qc,e(u) = q\ <7f, £ (w) = q = q' 



7 h+i t k y k+ \ fi 



70 =7 



Figure 5. The horizontal free path when ye T\q and y' e Ti,q \ 



and the contribution to G/ g(/l) of the interval [y, y'] is given by 

1 - W 1^/1, 

arctan tv - arctan y = • r- + O —z r- . (3.2) 

Y q(Kq + q>) 1 + y 2 \q 2 (Kq + q') 2 ) 



Since gg' > an d 2ye:F, e ^ ^ 1> ^ follows from ( 13.21 ) that the total contribution to G(j E (A) of intervals 
[y, y'] £ / with y < y' consecutive in Tq and y e T {t) is given by 



}fUA) = y y — L — +0 -\=s IQ (A) + o - 

CJ,eS > Zj Zj q(kq + q') 1 + y 2 \Qj WX ' \Q 



I.Q 

ft-<-lG<9A-+l 



When Y e similar bookkeeping with q' k = kq' + q, a' k = ka' + a, u k = provides the 



contribution 



G/ 3 ,' 0*) = y y (arctan y - arctan m a .) = Y Y — — — — 7 + ( — I . 

q' k <AQ<q' M q' k <AQ<q' M 

The situation is analogous to the one encountered in [5, Section 5]. Consider the "Farey triangle'' 
T = {(x,y) e (0, l] 2 : x + y > 1), the sets 



{(x,y) 



Q* = { (x, y) e R 2 



A-y 



x _ 

and for q e 24, fc > 1, the intervals 



> 0, 4 



A- 1 A- 1 



k k-l 



n[0,l), 



Denote 



/( 0) = (A - fe e ^, 1] = {e 9 ' : ( ei? , s 9 ') e Q ji _ 1 nf)c(l- £ g,l], 
= (1 - s^, A - ksq] = [eq : (sq, eq') e^nf)c(l - eq, 1]. 

, I — eq 1 

fk, q (q , a) = — — — — x , gi, 9 (jf,y) := A?(*> ? - W. « > °- 

q(kq + q') l + ^aY 

Assume first A > 2. Then min{/fc : Cl k n T} # 0} = [A] - 1 > 1, min{/t : 7 k £ 0} = [/I] > 2, and we 
have 

CO CO 

S/.eW) = Zj X fk,q(q'' a ) = Zj 2 ("S l,Q,k(^ <?) + Ti,Q, k (A, 



with 



Si,Q, k (A,q)= ^ f Kq (q,a), T mj! ,{A,q)= ^ fk-i, q (q' ,a). 

(a,q')eqIxQjl]l (cq'yeqlxQJ^ 
-aq'=\ (mod g), l\{q-a) -aq' = \ (mod g), £\{q-a) 
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Taking x = q' ', y — q - a — iz e q{\ - /), and I the multiplicative inverse of £(mod q), we can write 

Si,Q,k{A,q)= 2_j fk,q{x,q-y)= ^ gk, q (xjz), 

(x,y)eeQj'^xq( 1 -/), fe^g^ x (?/ f X 1 "O. 

.vv= I (mod g), C\y gcd(x,q)= 1, *zrf (mod g) 

Ti,g,k(A,q) = 2 g k -i, q (x,£z). 

(x,z)eQj'^x(g/e)(l-Il 
gcd(x,q)=l,xz=( (mod ij) 

Since (g, x) € QO/;, we have g - + 1)<? < < ^ - kq, so (A - l)Q ^ AQ - q < kq + x < t and one finds 

1 3 

qQ q 2 Q 



Since ^ | {q — a) and gcd(g, q - a) — 1 we have gcd(^, q) — 1 . The length of each of the intervals Qjf 1 
and Qj\ is less than g, so we can apply Lemma[3]with T = [Q Cl ] to find 

(i) 



Si,Q.k(A, q) = —5- ff gk,q(x, (z) dxdz + S ( , 

q 2 JjQjiyni-n 

* (3 3) 

dx |£ c t) 



with 



g/ <p(q) Q-q f AQ ~ kq dx 
£ q 2 Q J Q- q x + k< 

1 q^-^r- 

ZkU «« ^ q i+d + T^- qh s + Q \ rQ « Q^q-i + ° + ^ . (3.4) 
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Tl,Q,k(A 



c, y{q) Q-q f Q dx t 

( q 2 Q JAQ-k q x + {k-\)q \ 



A similar argument leads to 

ci tp(q) Q-q r Q 

lAQ-kq 

But 

r AQ - k " dx r Q dx r AQ du + r Q+ ^ du 

•J Q-q X + kq J AQ-kq X + (k - l)q jQ+(k-l)q u J AQ-q U 

hence ([3T3]»— (EIS) yield 



In 



1 I 
AQ-q' 



(3.5) 



S i,Q,k(A, q) + T IAk (A, q) = — 2 

I q 



ci <fi(q) Q-q 



In- 



AQ 



2c, -1 -i+i5 



+ O s AQ^~ L q 



Q- 



(3.6) 



Q *Q-q 

The intervals //. are disjoint, so when summing over k and q e Qlf, (or in a smaller range) we are actually 
summing over q e [ 1 , Q] . This way in d3.61 > the error will sum up to 



O, 



s.t 



G^-'jy^ + Q- c - c ' 2 -1 = o«(Q 5 " e(c,Cl) )> 



<]<Q 



while the main term will sum up to 



9<G q > 



MtjiQ) = y X — W(q), 



where 



gcd(f, 9 )=l 



wt ^ Q ~ 1 1 AQ 
W(q) = — In ■ 



qQ AQ-q 
with IIWHoc «: ^ and <K ^. Lemma Q] now provides 



9 €[1,6], 



M £>7 (0 = 



c,C{{) 



W(q)dq + 



AAnQ 
Q 



and we proved 

Proposition 3. For any 5 > 0, uniformly in A on compacts of [2, 00). 

,1 



An identical formula holds for Gy Ie {A). 



r 1 - u 

Jo « 



ln^ 



A — u 



du + O^s- 6 ^). 



(3.7) 




Figure 7. The set U~ ^ n 7" when 1 < A < 2 
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When 1 < A < 2 the contribution to S i,q(A) of y with q + q' < AQ is 

oo 



(q,q')£Q(n t nT) 



"Yj Yj ( fk,q(q',a)+ fk-i,g(q',a)j 

cdf^il ( fl '9' )e « /x e- / fc, ) (a,q')eqIxQjf> 

aq'=-\ (mod q),{\(q-a) aq'=—\ (mod ^(^-a) 

+ 2 2 /!*(«'• 

a-l)fl<?<!2 (a^eqlxQjf 
gcd(f,?)=l , , , . . '? . 
=-1 (mod g), 

and the contribution of y with q + q' > AQ is 

U-l)2«?<e (a,q')eqIX(AQ-q,Q] 
gcd(t,q)=\ aq'=-l (mod g),4(?-a) 

Using |L/yU < i and ||V/ t , 9 |U < on (Q/g U QJ$ Uq ) X ql if k > 1, ||/ , ? |U < and 
l|V/o ;? IL ^ (^-i)g 2 Q on 2*^1 « x 9 ^' anc ^ summm g as m me case ^ > 2 we also obtain 
Proposition 4. For every <5 > 0, uniformly in A on compacts of (1,2], 



,(2) n , _ ciC(€) 



f ' — In _i_ rf M + CV( £ -^>>). (3.8) 

Jo M A — u v 7 



An identical asymptotic formula holds for G\ r (A). 

As a result formula d3.7t will also hold when 1 < /I < 2. 

Consider finally the case < A < 1. When q' > AQ we have > AQ for all k > 0. So the 
contribution of each y e T^fg with > AQ is in this case 

1 1 J 1 
arctan y - arctan y = — ■ + O 



qq' 1 + y 2 \ q 2 q' 2 J ' 
summing up to 



K 9 <g (a,q')eqlXj A ^ 99 1+ (f) 
gcd(f, ? )=l god( ? ', 9 )=l, f|(?-a) 



a^'=-l (mod q) 

Using the same estimates as in case A > 2, d3.9t leads to 



f ~" ^ 77~ -du + O s ,{Q 6 - B ^) 

I Jo m max{l - m, A] v 7 



(3.10) 



Finally the contribution of each y e T\ C) q with </ < AQ is 



1 - eq 1 

arctan fo - arctan y = • r + , , , 

qq' I + y z \<? <? 



summing up to 



d iq w= y y 



a-,i)e<4«e (fl,?')e?/x(G-?.-ie] 

gcd(f,«)=l gcd( 9 ',«)=l,fl(g-a) 
ag'=-l (mod 



(3.11) 



f ' I — " i n _jL_ du + O sl (AQ- x+s + g^CA-i.-c-cd) . 
Ji_ A u 1 - u v 7 
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From ( 13.10b and ( 13.1 U we infer, uniformly in A on compacts of (0,1], 
Gf Is (A) = Ci, (A) + D W (A) 

^W/r'l, 1 ... ,\ . „ . ^ + max f2c ,-L- c - c ,A (3 ' 12) 



. |JT I In ^-L- du-Aj + Os,c(AQ- l+s + gffwawPci-l.-"*^ 



whence 



Proposition 5. For every (5 > 0, uniformly in A on compacts of (0, 1], 

An identical formula holds for G ( I (A). 

4. End of the proof of TheoremQ] 
From Corollary [1] and Propositions [3] El H] we gather 

G t ,,, e (A) = c&M + Os,e(s- s+e(c ' c,) ), 
with repartition function G( given by 



G C (A) 



1- ^+A{€)Hi(X) if^6(0,i], 
\--^ )+ A(€)Hi{X) ifAe[|,l], 
^H 3 (A) ifAe[l,oo), 



where 



Hi(A) = A- ((2) + (lni)ln(l -A) + 2 -In du + 2 \ In du, 

J A 11 I — U Jl-A 11 I — U 

r 1 1 i r A l l 

= -A - (ln/l)ln(l - /I) + -In du- -In du 

Ji_ A u 1 - u Jo u 1 - u 

= -A, 

, r 1 1 - h i 

i/ 2 (^) = A - f(2) + (In A) + 2 In du 

Jj M 1 — U 



= 3/l-2-^(2)-(ln/l) 2 + 2(l - A)lnf- - l) + 2 f i In — !— 

\/l / Ja 11 I - u 



1 1 - u . /I 



^(i) = In d M . 

Jo u A — u 



This establishes part (i) in TheoremQ] Using also 

f 1 1 - ull 1 \ 1 f 1 1 -m 1 / 1\ / 1 

ff'(A)= \du = — du = — + 1- - In 1-- 

3 Jo u \A A — uj A J A-u A \ A) \ A 

it follows that the density of the repartition function G( is 



A> 1, 



gt {A) = -G' e (A) 



^ r) +A(£) ifAe(0,i], 
^+A(0(-3 + |-2(i-l)ln(i-l)) if Ae[±,l), 
l^(i-(l-l)ln(l-i)) if A e (!,«»). 



Part (ii) of Theorem [TJ can now be deduced from part (i) by a standard approximation argument 
fl] lUHISl based on r f , £ (w) * for tanwe/ with / c [0, 1] interval of length |/| x £ e . We skip the 
detailed proof, which can be easily reconstructed from some of the arguments in the next section. 
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5. The conversion from a honeycomb to a square lattice with congruence constraints 
In the situation of the honeycomb it suffices to consider oj e [0, f ]. The linear transformation 

R 2 3 (x,y) (x',y') = ( x - — , — | e R 2 

\ V3 V3/ 

maps the first sextant r+ := {{q + |, ^y~) '■ q, a e Z, < a < q] of the grid of equilateral triangles of side 
1 onto the first quadrant of the square lattice Z 2 = {(q, a) : a,q e Z). Elements of the subset F+ ex c r + 
of vertices from the honeycomb grid map to integer lattice points (q, a) with q £ a (mod 3) (see Figure 




Figure 8. The free path length in the honeycomb and in the deformed honeycomb 




Figure 9. Change of scatterers under the linear transformation T 



T also maps circular scatterers (q + |, ^y^) + e(cos0, sinf?) centered at (xo,yo) = (q + | , ^y^) to 
ellipsoidal scatterers (q, a) + s(cos 8 - ^ji , centered at (jc^, y' Q ) = (q, a). Denote u>q - arctan ^ and 

ai' = arctan y. Denote also S B = {<5(cos |, sin |) : \S\ < s). 

In the honeycomb with scatterers {(xq, yo) +S e : (xq, yo) 6 r^ ex } of same length, the "local" repartition 
function 

\co e arctan / : q^ x (to) > - \ , / £ 0, - interval, 



(5.1) 



of the "horizontal" free path length 



<? £ hex (w) = inf \q e N : 3a e N, 3(x , yo) e 



r hex / g + Q/2 aV3/2 \ 
+ ' \ cos a> sin a> I 



e (xo,.yo) + S t 



turns out to be closely related with Gsj tS (X), being estimated through the same approximation procedure 
as for the later when I c [0, 1] is a short interval of length |/| x s°, < s < 1. Indeed, the equality 



T-\x ,y' + S) = \x + 



, y' + 6 (y' + 6)^ 
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shows that T maps the oblique scatterer (xo,yo) + S s from the honeycomb onto the vertical scatterer 
C^O'}^) + ^ e f rom me square lattice with x' Q = y' Q (mod 3), and the line through the origin with slope 
onto the line through the origin with slope <t>(-^), where <1> is the bijection 



O : [0, 1] 



V3 



:V3 



2x 



2+fx' ' V3-x 

In the process the condition q^ ex (co) > ^ above is being replaced by q^^co') > - in the square lattice. 
The only difference arises from replacing expressions (with x' = q, y' Q = a and < 5 < -) 



arctan 



?0 + (5 



arctan — - — 

x x {) 1 + 



1 



q 2 + a 2 



that collect the contribution of angles to for G3j >e (A), by 



arctan <t> 



arctan Ol I = arctan 



(y' Q + <5)V3 



arctan 



VI 



2x' Q + y' Q + 6 2x' + y' Q 



2x' Q +y' +S 2x' n +y' n /. i y Yi y 
L+ \ 2 *'o+y'J 

5 \\ 1 



/ V3 \ 2 



1 +01 



\2x' +y' )) l + (y^Lf 



+ o U 



2x' 6 V3 
~ (2x^+^) 2 

_ 6q^j3 

2(q 2 + aq + a 2 ) \q 2 / 

The effect will only be on the main term, where cj — J } j^p, J £ [0,1], will be replaced by cj ex := 
Jo-'(/) 2(i+fL } ) I - "Tf J' ^ n wa y we obtain, uniformly in A on compacts in R + \ {1}, 

G£f (A) = cY x G 3 (A) + O s {s- d+e(c ' Ci) ). (5.2) 
The change of variable x = <£>(//) gives 



J [tano»o,tano)i] 



V3 r 

2 Jo 



an o»i ) 



dx 



'(tan Wo) 



X + X + 



- P 

1 Jtan 



d\JL 
fi 2 + l 



= 0)1- co . 



(5.3) 



In this case q^ ex (oj) and the free path length (co) are related (by the rule of Sines) by 

~ hex/ 



T^V) 



Sin ^ ^ex r/ , V3 q^io) 
q B (co) 



sin(? - co) 



so that 



'^(A) := 



w e arctan / : f^icS) > - 
e 



2 cos (| + w) ' 



r , 2/1 cos (5 + co) 

\co e arctan/ : q^ x (co) > - h 



:V3 



(5.4) 



Fix a»/ e arctan/. Using | cos (| + co) - cos (| + < \io — a>i\ <K e c , (12.2b . and c^ ex « e c , equalities 
Ol and (El yield 



P^(A) = c* ex G 3 



2/1 cos (| + co;) 



+ 6 (s~ i c 



2c , -S+9'c,ci) 



(5.5) 



Let e± = e ± (fa;o,e) as in Figure [9] Consider the case of scatterers (xo,yo) + S WOj£ , (xo,yo) £ F 
where S WOj£ denotes the segment {<5(cos |,sin|) : -e_ < £ < £+}. Let f £ ex (w) denote the free path 



hex 

+ 
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length and P^(X) = \{co e arctan/ : fg ex (w) > Since |cos a> ± - cos <sc e c and cos(£ + a>j) > \ 
there exists C\ > such that 

' -Ci£ f | <e_ <e + <e" :=e( -+C I e e |. (5.6) 



^COS(|+W/) / \COS(|+(jJ/) 

The inequalities f^„(w) < fg(w) < f^,(aj) and (15. 6t combined with formula d5.5t and j2.2\ lead to 

P ff(A) < 



Jai e arctan/ : r]f x (<j/) > - 



to € arctan / : f^, ex (w) > 



£'(cos(|+w / )-C 1 s c ) 

/ J \ 

hex 



- c 7 G3 



COS(| + W/) - C\S° 

2A cos(|+W/) 



-J + O g (e 2c + s- 6 ^^) 



V3 cos(| +w / )-C 1 e e 

and to a similar lower bound for P 7 h ™(/1), and so we get 

p*f(A) = cJ ex G 3 |-^J + (5 (£ 2c + g-^^.'')). (5.7) 



The trivial inequality \t^(oj) - fj((t))\ < 2e and ( 15.71 ) now provide the formula 

ph-^) = C J-G 3 |^|J + O s (s 2c + s - 6+B ^ (5.8) 

for the repartition function P / h ™(/i) = \{a> e arctan/ : t^ x (co) > |}| of £T £ hex . 

Finally we choose a partition [0, = uf =1 Ij with intervals 7 ; of equal size ^ - e c . Applying J5.8t 
to each individual interval Ij with c = Ci = 1 and summing over 7 we find 

N 



2) hex 
' [0,1/' 



- *.» - g «r«> (^) + - H^) + °"( sH) ' <59) 

Theorem|2]now follows immediately from (15. 9t and obvious symmetry properties of the honeycomb. 
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